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We use random-matrix theory and supersymmetry techniques to work out the two- 
point correlation function between states in a hierarchical model which employs 
Feshbach's chaining hypothesis: Classes of many-body states are introduced. Only 



states within the same or neighboring classes are coupled. We assume that the 
density of states per class grows monotonically with class index. The problem is 
mapped onto a one-dimensional non-linear sigma model. In the limit of a large 
number of states in each class we derive the critical exponent for the growth of 
the level density with class index for which derealization sets in. From a realistic 
modelling of the class-dependence of the level density, we conclude that the model 
q ■ does not predict Fock-space localization in nuclei. 
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Statistical methods play an important role in the quantal many-body problem. 
We recall the importance of the Gaussian Orthogonal Ensemble of random 
matrices (the GOE) for the description of fluctuation properties of the energies 
and wave functions of neutron resonances and other data on nuclear energy 
levels located a few MeV above the ground state. The use of the GOE is not 
restricted to nuclei, of course, and this ensemble is equally important in other 
many-body systems like atoms or quantum dots [1]. 

It is well known that the GOE does not furnish a completely adequate de- 
scription of stochasticity in any of these systems. This is because the forces 
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between the constituents (nucleons in the case of nuclei, electrons in the case 
of atoms and quantum dots) are predominantly of two-body type. Then it is 
more natural to assume that the matrix elements of the two-body interaction 
are random variables. Both in nuclear and in atomic physics, there is empiri- 
cal evidence that this assumption is viable [2,3]. The random-matrix ensemble 
which implements this assumption in large shell-model spaces (the two-body 
random ensemble) differs from the GOE: The number of independent two- 
body matrix elements is generically much smaller than the dimension of a 
typical shell-model matrix, while for the GOE, the number of independent 
matrix elements is proportional to the square of the dimension of the matrix. 
Numerical studies have indicated long ago that this characteristic difference 
is immaterial for spectral fluctuation properties [4,5]. However, these stud- 
ies were limited to small matrix dimensions. Therefore, the question persists 
whether the GOE and the random two-body ensemble are fully equivalent. 

The question has resurfaced with the recent observation that in quantum 
dots, a statistical model for the two-body interaction allows, in principle, for 
localization in Fock space [6]. Numerical studies seem to support this proposal, 
at least at sufficiently high energies [7], while near the Fermi energy no evidence 
for localization was found [8]. It is then only natural to ask whether a random 
two-body interaction in nuclei might likewise lead to Fock-space localization 
in nuclei. Such localization would have exciting implications for level statistics 
and wave function properties. Moreover, Fock-space localization could lead to 
novel aspects of nuclear reaction dynamics. 

Localization is a wave phenomenon discovered by Anderson [9] . It occurs when 
non-interacting electrons move diffusively under the influence of impurity scat- 
tering through a conductor. The Anderson tight-binding model describes this 
process in terms of a lattice with random on-site energies distributed uni- 
formly over an interval of length 2W and with fixed hopping matrix elements 
V connecting neighboring sites. Localized eigenfunctions of the tight-binding 
Hamiltonian are not spread uniformly over the entire lattice but possess an 
envelope which decays exponentially at large distance. Anderson localization 
occurs not only in the passage of Schrodinger waves through a disordered 
medium but likewise in the passage of light through a medium with a disor- 
dered index of refraction [10]. 

The concept of localization can be extended to many-body systems. This was 
first suggested in quantum chemistry [11]. Logan and Wolynes observed a 
close topological similarity between a statistical description of the many-body 
problem in polyatomic molecules and the tight-binding model. An analogous 
similarity exists in quantum dots [6]. In the present paper, we point out that a 
corresponding similarity also exists in a Fock-space description of many-body 
systems like nuclei or atoms with random two-body interactions. We investi- 
gate the question whether in these systems, localization is expected to occur. 
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The occurrence of localization can either be established by very extensive 
numerical calculations or analytically. Here, we choose the second approach. 
First, we must seek a model which allows us to pose quantitative questions. 
Optimally, we would want to consider a Hamiltonian which is the sum of one- 
body and two-body operators. The one-body terms would describe a nearly 
degenerate set of single-particle or shell-model orbitals, and the two-body 
terms would contain the random two-body matrix elements. This kind of 
model does not allow for an analytical treatment, however, and it is difficult 
to see how numerical calculations on the required scale could be performed. 
Our work is, instead, based on a Fock-space model commonly used in the sta- 
tistical theory of nuclear reactions. We recall that in nucleon-induced nuclear 
reactions, precompound reactions are important whenever the equilibration 
time of the compound system is comparable to or even larger than the decay 
time of the system. This situation occurs typically at bombarding energies of 
several 10 MeV. Then, the standard compound-nucleus model (which uses the 
GOE) fails to apply and is replaced by a hierarchical statistical model: The 
incident nucleon creates a sequence of m-particle (m — l)-hole states where 
the integer m increases by one unit in each two-body collision of the projec- 
tile with one of the target nucleons. The density of m-particle (m — l)-hole 
states at fixed excitation energy grows very strongly with increasing m, see 
Eq. (27) and the text following it. Therefore, a dynamical description of the 
process is out of the question and statistical concepts are used instead. We 
define classes of states. The states in class m are composed of the m-particle 
(m — l)-hole states. Within each class, the two-body interaction mixes the 
m-particle (m — l)-hole states. We assume that the resulting matrix problem 
in each class m can be replaced by a GOE. The GOE's pertaining to different 
classes (different m values) are uncorrelated. For each m, the one parameter of 
the GOE is fixed by the density of states in that class. Classes differing in m 
by one unit are connected by the elements of the two-body interaction which 
are taken to be independent Gaussian distributed random variables. This last 
assumption corresponds to Feshbach's "chaining hypothesis" . It embodies the 
two-body character of the interaction and defines a hierarchical model. Within 
this model, we ask whether the eigenstates of the full Hamiltonian are local- 
ized. We now turn to the details of this model. We closely follow the work of 
Ref. [12]. 

We introduce a basis |m/z) of states where m is the class index introduced 
above, and where /z with fi = 1, . . . , N m is a running index. The Hamiltonian 
reads 



H = \mfi)H m ^ nu (nu\ (1) 

mnfii/ 
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and 



is real and symmetric under exchange of (m/i) and (W). We assume that all 
matrix elements are Gaussian distributed random variables with the following 
first and second moments, 



H m /j,,nu $rnn$iiv h m , (3) 



<5 ' H m ^ nv 8 H-m' fj,' ,n' v' -^mn mm 1 ^ nn' ^ jipJ ^ vv' $mn' $nm' & \w' &vjj/) ■ (4) 

The bar indicates ensemble averaging. For fixed m, the matrix H m ^ mv belongs 
to the GOE. The matrices describing different classes m are seen to be mi- 
correlated. In order to keep the GOE spectrum finite in the limit N m — > oo, 
the diagonal terms M mm must scale as iV" 1 . We consider the weak coupling 
case [12] where the non-diagonal elements M mn (m 7^ n) are suppressed by 
an additional factor N^ 1 , 

M mn = §^5 mn + - 5 mn ) . (5) 



Here \ rn and \ mn are the strength of the matrix elements in class m and 
the coupling strength of states in different classes, respectively. In keeping 
with the remarks made above, we assume that \ mn vanishes for \m — n\ > 2. 
The spectrum of the N m states in class m has the shape of a semicircle with 
center h m and radius 2A m . We proceed by assuming that all semicircles have 
identical centers, h m — for all m, and identical radii, \ m = A for all m. The 
remaining parameters N m allow us to account for the fact that the local level 
densities p m (E) at energy E differ in the classes. We work at the centers of 
the semicircles and put E = without loss of generality. Later we need the 
inverse g of M given by 

9mn = ^S mn - (1 - 5 mn ) + O (^^) • (6) 

With D(E ± ) = E — H ± irj and the ensemble-averaged two-point correlation 
function /C defined by (m 7^ n) 



lC m „,nv = (m/i| J D- 1 (E+)|n^)(n^| J D- 1 (^-)|m/x) , (7) 



the quantity of central interest is the class average of IC, 

K-mn — "77 TJ~ ^ ^mfi-nu ■ (8) 
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If K, mn vanishes exponentially with increasing distance \m — n\, the eigenstates 
of H are (exponentially) localized. 

The explicit dependence of JC mn on \m — n\ can be found with the help of the 
super symmetric non-linear sigma model [13-15] the use of which is by now 
completely standard. Moreover, a derivation taylored to the present problem 
may be found in Ref. [12]. Therefore, we restrict ourselves to the essential 
steps and emphasize those aspects which are specific to the present problem. 
We use the notation of Ref. [15]. The correlation function is expressed in terms 
of the generating function Z, 



/V/rn.r). 



(9) 



where J^ 5 denotes the (1,5) element (with respect to the supersymmetry in- 
dices) of the auxiliary J-field. Similar notation is used in Eq. (15). Here 
we have averaged the generating function over the Gaussian distribution of 
H mix nil . After eliminating the quartic dependence on the original integration 
variables by means of a Hubbard-Stratonovitch transformation, we find for Z 
the expression 

Z(E,J) = 

J d[a] exp ^~ E A ' 9mn trg(a m a n ) - X - trg m /i Q In M{ J^j . (10) 

The fields a m are 8x8 graded matrices and d[a] = Y[ m d[o'm\. Moreover, 

Af(J) = E + ir) + J-X, E = {A5 mn ^<' ? }. (11) 



The integral in Eq. (10) is evaluated by means of a saddle-point approximation 
by varying the <7 m 's. In view of the smallness of the terms with m ^ n in 
Eq. (6), we obtain a separate saddle-point equation a m (E — Xa m ) = A for each 
class to. The solution yields the semicircle for the level density. Integrating over 
the massive modes in the limit where N m — > oo leaves us only the Goldstone 
modes <r^. We focus attention on the correlation function fC(m) = K\ m which 
takes the form 

1 r M /M-l \2 \ 

K{m) = 2^J n^HK G ) 51 (^) 51 ex P E ^trg(afaf +1 ) . (12) 

k=l \ i=l I 



Here M is the total number of classes. In the case of weak coupling we have [12] 

4% = 2tt Pm \TJ 2tt Pn (13) 
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where we wrote M mn = V mn 2 and where p m = N m /ir\ is the density of states 
in class m at energy E — 0. The form of the expression in the exponent 
of Eq. (12) displays the chaining hypothesis (only states belonging to next- 
neighboring classes are coupled by non-vanishing matrix elements). 

To proceed, we omit from now on the index G on the sigma fields and follow 
Refs. [16,17]. The aim consists in replacing the summation over classes by an 
integration (continuum limit). We put 



X 2 

A u+i 



1 



4A 2 4e 



a 



i,i+l 



(14) 



with e = 1/C and both ( and a iji+ i ^> 1. This is justified because the left-hand 
side of Eq. (14) is much larger than unity, see Eqs. (26) and (27). We write 
the correlation function in the form 

JC(m) = ^2 / d°id(Tm erf W((T X , a m ; 1, m) Y(a m ; m) (15) 



where 



ra-l 



W(a m ,a n ;m,n) = J J| da k exp 

k=m+l 
r M 

Y(a m ;m) = J J[ da k exp 



k=m+l 

Obviously, W obeys the equation 



n-1 



J2 a v+i tr g( ff » ff i+i) 



i=m 
M-l 



-£ tig (aia i+1 ) 



(16) 



W(a m ,a n+1 ;m,n + 1) 

= J da n W(a m ,a n ;m,n)exp 

We use the smallness of e to take the continuum limit and replace the sum- 
mation over classes by an integration over the continuous variable t = m/(. 
Then, 



—a„,„+i trg(o- n cr n+ i) . (17) 



W(<7i,<7 m ; l,m) 
Y(a m ;m) 



a(t) 

Y(a(t);t) 



and Eq. (17) becomes an integral equation. This integral equation is equivalent 
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to a modified heat equation, 



a(t) d t W(a, a; t) = A a W(a', a; t) (18) 
with the initial condition 

lim W(a', a; t) = 5(a' - a) . (19) 
Here A a = d 2 /da 2 [16,17]. 

In previous work [13,17], the prefactor in the exponent of Eq. (12) was a con- 
stant (independent of the class index i), and the solution of the heat equation 
did display an exponential decay signalling localization. In the present case, 
a(t) does depend upon t because the level densities p m depend strongly upon 
m. This fact raises the following questions: (i) Is localization affected or even 
destroyed by a monotonic increase of p m with m and, thus, of a(t) with t? (ii) 
If so, can we determine the critical exponent f3 in a(t) oc t 13 which marks the 
transition from the localized to the delocalized regime? (iii) What are the im- 
plications of our findings for localization in nuclei, i.e., for a realistic modelling 
of the t-dependence of a(t)l 

To answer these questions, we observe that we retrieve the standard heat 
equation by a change of variable, 

a(t)d t = d s . (20) 

In terms of the rescaled function 

W(a u a(s);s) = Wfa, a(s(t)); s(t)) = W(a u a(t);t) (21) 

this leads to 

d s W(a',a;s) = A a W(a',a;s) . (22) 

We choose s(t) such that s(0) = which ensures that W also obeys Eq. (19). 
The transformation (20) reduces our problem to the heat equation solved in 
Ref. [13]. The rescaled propagator 

1C(s) = ^2 / da 1 da( S )al 1 W(a 1 ,a( S y,s)a 51 ( S )Y(a( S y,s) (23) 
can be worked out by introducing the Fourier transform K{k) of Eq. (23). 
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This yields 



K{8) = 



7f 

16CA 2 



:?)*'•-(?: 



(24) 



which is valid for s ^> (. The function /C(t) is obtained by solving for s(t) the 
differential equation 

df(0 _ , . 

~a(f)' 1 J 



and replacing s by s(t) in Eq. (24). For a power-law dependence, a(t) oc t' 3 , 
Eq. (8) yields s oc t 1_/3 . Combining this with Eq. (24), we see that f3 — 1 is 
the critical exponent which marks the transition from localization to dereal- 
ization. The existence of such an exponent is intuitively clear: A monotonic 
increase of a(t) with t signals an ever increasing phase space. The resulting 
drag on the system overcompensates, for (5 > 1, the tendency of the system 
to form localized states. 

In order to apply this result to the case of nuclei, we must determine the 
dependence of a(t) on t. We recall that 

a(t) <^> a m ^ m+ i = (tt p m V m 

,m+l Pm+1 )/C = (vr 2 



where r^ m+1 is the spreading width for transitions from class m to class 
m + 1. The width r^. t ^ m+1 approximately scales with 1/m [18]. The density 
of states in class m at excitation energy E is [19] 



Pm(jE) "m!(m + l)!(2m)!AUJ (2?) 



where A is the mean level spacing of the single-particle states. For fixed E, 
p m {E) grows strongly with m up to a maximal value of m and then drops 
quickly. In the present context, we are not interested in this cut-off because it 
might mask the existence of localization. For this reason, we take the limit of 
high excitation energy E/A 3> m. Then, we can neglect the factorial factors in 
Eq. (27) as well as the 1/m dependence of Y m ^ m+ i and obtain approximately 
a m , m+1 oc (E/A) 2m . This yields 

s(t) & s m oc 1 - (E/A)- 2m . (28) 



We observe that a mtm+ i grows with m much more strongly than linearly. As 
a consequence, s m is bounded from above and approaches asymptotically a 
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finite constant. Thus, the correlation function JC(t) tends to a finite value as 
t — > oo. We conclude that the interaction spreads the states in class 1 over 
the entire Hilbert space. 

The results of the model investigated in the present paper negate the possibil- 
ity of localization in a many-body Fermi system like the nucleus. This finding 
must be contrasted with the results of Ref. [6] where the chain of m-particle 
(m — l)-hole states was mapped onto a Bethe lattice, and localization was 
predicted. Such mapping completely neglects the interaction between states 
within each class m. In contradistinction, our use of a GOE for the states 
within each class may overemphasize the role of the interaction between states 
within each class m, although we believe the present model to be closer to re- 
ality in nuclei than the model using a Bethe lattice. As mentioned above, in 
the numerical work of Leyronas et al. [7] evidence for localization was pre- 
sented. In that work, all those m-particle (m — l)-hole states were taken into 
account the unperturbed energies of which lie in some small neighborhood of 
the energy E considered. It is not clear to us whether this constraint might 
not affect the localization properties of the system. 

In summary, we have shown that in the framework of our hierarchical model 
localization is destroyed whenever the density of states increases more strongly 
than linearly with increasing complexity m of states. In the case of nuclei, the 
model does not predict localization. 

T.R. thanks C. Mejia-Monasterio for beneficial discussions at the early stages 
of this work. 
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